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TMultidimenaloDal  scaling  procedures  provide  a new  means  of  studying  spatial  organi- 
zation of  brain  electrical  activity.  A three-staged  analysis  process  Is  proposed  and 
tested  on  human  electroencephalographlc  dataa  Simple,  high  accuracy  solutions  with 
anatomically  meaningful  dimensions  were  achieved.  The  procedure  Is  not  limited  to 
the  application  presented,  but  is  of  general  utility. 
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Multidimensional  scaling  prox.  res  provide  a nsw  means  of 
st'tdying  spatial  organi7ation  of  brain  electrical  activity.  A three-staged 
analysis  process  is  proposed  and  tested  on  hum.?j\  eleotroenccphalographic 
data.  Simple,  high  accuracy  solutions  with  anatomically  meaningful  dimen- 
sions were  achieved.  The  procedure  is  not  limited  to  the  application  presen- 
ted, hut  is  of  general  utility. 


INTRODUCTION 

In  tite  elec^rcphyeiologlcal  analysis  of  brain  fimction,  elegant  and 
powerful  techniques  have  been  developed  to  study  the  temporal  pattern  of 
activity  at  a sii:glo  recording  site,  but  rather  less  progress  has  been  made 
in  solving  the  more  difficult  problems  of  anal3'7ing  the  spatial  patterning  of 
activity  at  a number  of  simultaneously  recorded  sites  (1).  Ultimately  much 
of  our  imderstandirg  of  the  organiration  cf  brain  functions  must  depend  upon 
such  configurational  analyses  (2).  In  the  present  paper,  I report  that  re- 
cently developed  techniques  of  multidimensional  scaling  provide  useful  an- 
alitic  toe's  for  the  extraction  of  information  about  the  fuactlonal  relations 
between  the  sources  of  simultaneously  recorded  multichannel  electrophys- 
iological  data,  using  the  multichannel  electroencephalogram  (££G)  as  an 
example. 

The  procedure  involved  consislo  of  three  logically  distinct  steps. 

In  the  first,  multi-channel  elcctrophysiological  data  is  acquired  from  the 
preparation,  inspected  for  artifacts,  preprocessed  if  necessary,  and  then 
stored  for  analyeis.  In  the  second  step,  some  aspect  of  that  data  is  selec- 
ted by  appropriate  featimo-selectlon  and  proximity  estimation  algorithms 
which  are  applied  to  all  chancels  of  data  in  a pairwise  fashion.  The  product 
of  this  step  is  a numerical  estimate  of  the  similarity  or  dissimilarity  of 
each  pair  of  channels  v/ith  respect  to  the  features  of  interest.  In  the  third 
steo,  the  rosuiting  previTnlty  matrix  is  analysed  by  a multirlimcnsional  scal- 
ing procedure  to  dntormiiic  both  the  number  of  d;‘m'?r.oicns  necessary  to  re- 
produce the  measured  oimilarlty  and  the  configuraticn  of  the  channels  in 
that  multidimensional  space. 


Nonmctric  multidJmeoaional  scaling  analysis  attempts  to  derive 
the  configuration  of  pointa  In  on  N-dimensiccol  Euclidean  space  such  that 
the  diotances  between  points  Ln  the  configuration  are  as  nearly  as  possible 
a monotonic  function  of  the  data.  As  an  iUimtrative  example,  a proper  map 
of  major  American  cities  may  be  reconatructod  from  a half  matrix  of  air- 


>Vhife  Section 
Tjtt  Section 


□ D 


line  distances  between  these  cities  when  solved  in  two  dimensions  (3).  In 
the  present  investigation,  the  object  Is  not  to  reconstruct  the  spatial  config- 
uration of  electrodes  placed  upon  the  subjects'  heads,  but  to  use  the  scaling 
method  to  determine  the  functional  configuration  of  the  brain  regions  from 
which  the  EEC  was  recorded.  Function  proximity  reflects  similarity  in  elec- 
trophysiological  activity. 

METHOD 


EEC  was  recorded  from  10  normal  young  adults  who  sat  with  their 
eyes  closed  in  a quiet  room.  Monopolar  recordings  were  taken  from  both 
left  and  right  cer^ral  hemispheres  in  the  frontal  (F3  and  F4),  central  (C3 
and  C4),  parietal  (P3  and  P4)  and  occipital  regions  (01  and  02),  with  the 
linked  ears  as  reference.  The  letter-digit  pairs  refer  to  standardized  sites 
for  EEC  recording,  which  are  illustrated  in  Figure  1 (4).  EEC  was  ampli- 
fied by  a Beckman  R-411  polygraph  (bandpass:  0.16-30  Hz,  gain:  SO  micro- 
volts/cm). Five-sec  samples  ci  EEC  were  digitized  at  10  msec  intervals 


and  stored  on  the  disc  memory  of  the  laboratory  computer  system  (Hewlett- 


01  02 


Figure  1.  Electroencephalograi^c  recording  sites  as  designated  by  the 
conventions  of  the  International  10-20  System  (4). 


Packard  2116B/2879A).  As  each  S-sec  sample  was  acquired,  the  paper  re- 
cord from  the  polygraph  was  examined  for  eye-movement  and  muscle  arti- 
facts. All  samples  containing  visible  artifacts  were  rejected.  l>velve  arti- 
fact-free samples  were  obtained  for  each  subject.  Figure  2 presents  a typi- 
cal 8-channel,  5-sec  digitized  EEG  record. 
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For  each  subject,  the  functional  prcxin'itj'  of  each  pair  of  EEG 
cbannelc  was  estin  atec'  with  respect  to  sharcc'  synchronous  waveforms  usin^ 
the  polarity  coincidence  corx'elatitm  (PCC)  coefficient,  which  is  computed  by 
countinc  the  total  number  of  instances  In  which  the  polarities  of  the  two  time 
series  are  of  the  same  sign,  subtracting  from  that  count  the  number  of  occa- 
sions on  which  the  polarities  are  of  opposite  sign,  and  dividing  this  difference 
by  the  total  number  of  pairs  counted.  This  statistic  ranges  between  +1.0  for 
perfectly  synchronous  signals,  and  -1. 0 in  the  case  of  perfect  phase  rever- 
sal, with  0. 0 indicating  a lack  of  association  between  the  two  time  series. 

The  PCC  coefficient  is  trigonom ically  related  to  Pearson's  product-moment 
correlation  coefficient  for  Gaussian,  ergodic  randomi  signals  (5),  and  also, 
therefore,  varies  with  the  value  of  the  conventional  crona  -cori-olation  fimc  - 
tion  computed  at  ? ero  time  delay. 

Cince  negative  values  of  the  FCC  seldom  are  observed  in  elec- 
troencephalographic  recording  (6),  FCC  may  be  considered  to  reflect  the 
relative  amount  of  synchronous  activity  between  a pair  of  EEG  channels. 

TWO  EEG  channels  that  exhibit  highly  synchronous  activity  are  considered  to 
be  functionally  close  to  each  other  whereas  a pair  of  channels  that  show  lit- 
tle synchrony  are  considered  to  be  functionally  distant. 

Harmony  and  her  coworkers  (6)  have  shown  that  the  FCC 
Is  a stable  paran  eter  of  the  multichannel  electroenceplialogram,  varying 
little  within  a particular  individual  ovei:  repeated  tests.  Further,  the  value 
of  the  PCC  is  sensitive  to  at  least  some  changes  in  the  functional  state  of  the 
brain;  it  typically  increases  during  photic  driving,  for  example  (6).  More- 
over, abnormal  values  of  PCC  have  been  reported  in  children  with  learning 
disabilities,  suggesting  that  the  measure  partially  reflects  aspects  of  brain 
function  related  to  cognitive  processing  (7).  This  idea  receives  som^e  addi- 
tional support  in  that  Callaway  and  Karris  (8),  using  a similar  measure  of 
cortical  coupling,  have  shown  hen.ispbere-speciffc  changes  occurring  in  ver- 
bal and  spatial  inform  ation-processing  tasks. 

Applying  the  FCC  statistic  to  all  pairs  of  FEG  channels  resulted 
in  an  O-by-0  halfmatrix  of  PCC  values  for  each  subject.  These  proximity 
matrices  then  were  analysed  for  dimensionality  and  for  the  configuration  of 
channels  within  that  dimensionality  using  a general  !?-way  program  for  mul- 
tidimensional scaling,  KYST  (£).  ntandard  options  were  employed. 


in  Fi[:ure  3.  In  these  solutions,  points  that  cluster  together  are  functionally 
similar,  whereas  points  separated  by  greater  distances  in  the  solutions  are 
functionally  more  dlssiirilar.  Froir  these  configurations  it  appears  that  or- 
derly information  concerning  functional  proximity  of  brain  regions  is  con- 
tained within  the  multichannel  electroencephalogram. 

Although  the  configurational  solution  may  be  freely  rotated  in  2- 
way  Euclidean  scaling  procedures  such  os  KYCT,  in  the  present  example  the 
rotation  to  principal  con  ponents  achieved  automatically  by  the  program  was 
sufficient.  The  two  dinr  ensions  that  appeared  are  anatomically  meaningful 
and  may  be  labelled  anterior-posterior  and  left-right  respectively.  On  these 
(hmensions,  with  the  exception  of  a reversal  in  the  occipital  channels  of  Tub- 
Ject  6,  every  channel  is  appropriateiy  placed  with  respect  to  the  other  chan- 
nels given  the  anatomical  arrangements  from  which  the  electroencephalogram 
v/as  recorded  (10).  Put  the  distances  in  these  configurations  do  not  represent 
anatomical  distances  which  may  be  measured  in  millimeters,  but  rather  re- 
flect functional  distances,  which  are  measured  in  the  units  required  by  the  al- 
gorithm generating  the  original  proximity  matrices.  In  this  case  the  unit  of 
distance  reflects  the  degree  of  synchronous  EEC  activity  at  the  various  sites 
measured. 


These  ^-dimensional  configurations  represent  quite  accurate  fits 
to  the  data.  In  the  KYCT  procedure  for  multidimensional  scaling,  the  config- 
urations are  altered  to  minimize  "stress",  a measure  of  the  badness  of  the 
fit.  Ctress  as  used  here  is  defined  as  the  ratio  of  the  root  mean  squared  er- 
ror in  estimating  the  original  proximity  data  from  the  computed  dimensional 
configuration  of  distances  and  the  root  mean  square  of  the  distances  in  the 
configuration.  Truska.l  (11)  regards  a final  stress  value  of  .05  as  represent- 
ing a "good"  and  .035  as  an  "excellent"  fit  to  the  data  by  the  computed  ccm- 
figuration  in  many  situations.  By  this  standard,  the  fit  of  the  configurations 
shown  in  Figure  3 is  exceptionally  accurate,  the  stress  values  of  these  con- 
f4:.urations  ranged  between  .004  and.  .010  for  individual  subjects. 

V ith  respect  to  the  dimensionality  of  the  solution,  the  2-dimen- 
sional solutions  presented  appear  to  be  most  adequate.  The  stress  values 
associated  with  1 -dimensional  solutions  were  about  twenty  times  as  great  as 
those  for  the  2-diir«itlonal  solution,  raiding  between  .080  anc'  . 156  for  indi- 
vidual subjects.  V’lth  only  8 points  in  a configuration,  a 3-dimensional  solu- 
tion is  not  warranted  (C). 

An  examination  of  the  plot  of  computed  distances  in  the  configura- 
tion against  the  corresponding  proximities  in  the  PCC  matrix  (Shepard  dia- 
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grair)  reveals  a linear  or  near-linear  relation  for  every  subject.  This  in- 
dicates that  a metric  analysis  of  these  data  could  be  employed. 

DISCUSSION 

The  importance  of  these  findings  is  as  follows:  First,  multidi- 
mensional scaling  procedures  provide  a method  for  clearly  visuali7ing  and 
comprehending  the  functional  organization  of  electrophysiological  activity 
as  it  simultaneously  appears  at  a number  of  recording  sites.  This  provides 
a solution  to  one  of  the  most  difficult  problems  in  multichannel  electrophys- 
iological analysis,  representation  of  the  data  in  a eoneeptnally  rno.-iningfiil 
form. 

Second,  this  method  is  not  limited  to  the  study  of  polarity  coin- 
cidence of  the  human  EEC.  Any  source  of  multichannel  electrophv Biologi- 
cal data  may  be  employed,  providing  that  a substantively  meaningful  proce- 
dure for  estimating  functional  proximity  may  be  found.  One  obvious  appli- 
cation would  be  in  the  study  of  multiple-channel  averaged  evoked  potentials, 
where  a correlation  statistic  could  be  used  to  assess  waveform  similarity 
of  the  potentials  and  multidimensional  scaling  to  provide  a map  of  function- 
ally related  recording  sites.  Similar  approaches  could  be  applied  to  the 
study  of  amplitude  variations  in  extracellular  recordings  of  multiple  unit 
activity  in  any  brain  region.  Coherence  measures  of  cross-spectral  analy- 
sis (12)  also  provide  a proximity  measure  suitable  for  processing  by  multi- 
dimensional scaling  techniques. 

Third,  differences  between  configurations  may  be  explored  us- 
ing 2-way  multidimensional  scaling  programs  such  as  INCSCAL  (13).  This 
procedure  requires  linearity  in  the  data.  INBSCAL  operates  on  a set  of  re- 
lated proximity  matrices  and  extracts  a common  configuration  for  the 
group.  Individual  configurations  are  achieved  by  altering  the  relative  di- 
mension weights  of  the  group  solution  to  fit  individual  matrices.  This  pro- 
cedure is  useful  both  for  assessing  individual  differences,  as  in  evaluating 
multichannel  EEC  in  normal  and  neurologically  disordered  populations,  and 
for  studying  the  effects  of  various  experimental  treatments  within  individ- 
uals. 

Fourth,  the  substantive  results  of  the  present  experiment  indi- 
cate the  presence  of  a strong  pattern  of  two  dimensional  spatial  organiza- 
tion in  human  electroencephalogram  under  resting  conditions.  The  low 
stress  values  associated  with  the  solutions  obtained  for  each  of  the  subjects 
indicates  that  the  configurations  obtained  in  two  dimensions  fit  the  actual 


dat^  quite  closely.  The  fact  that  the  principal  components  rotation  of  the 
final  two  dimensional  configurations  resulted  in  an  anterior-posterior  and  a 
left-right  axis  for  each  of  the  10  subjects  suggests  that  these  solutions  re- 
flect natural  properties  of  the  data. 

For  these  reasons,  it  appears  that  multidimensional  scaling  pro- 
cedures may  provide  one  solution  to  the  neurc^hysiological  problem  of  sen- 
sibly analy?ing  and  comprehensibly  organizing  multiple- channel  electro- 
physiological  data. 
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